AD-A147  501  AN  ANALVTICAL  MODEL  OF  PERIODIC  HAVES  IN  SHALLOW  WATER  1/1 
<U)  AERONAUTICAL  RESEARCH  ASSOCIATES  OF  PRINCETON  INC 
NJ  H  SEGUR  ET  AL.  JUL  84  ARAP-TM-84-12 
UNCLASSIFIED  N00014-80-C-0048  F/G  12/1  NL 


OUC  FILE  COPY  AD- A 147  501 


A.R.A.P.  Tech.  Memo.  8*4-12 


AN  ANALYTICAL  MODEL  OF 


PERIODIC  WAVES  IN  SHALLOW  WATER  J 

0,Q strict  - 

by 


Harvey  Segur 

J 

Aeronautical  Research  Associates  of  Princeton,  Inc. 

50  Washington  Road,  P.0.  Box  2229 
Princeton,  NJ  085*10 


and 


Allan  Finkel 


Thomas  Watson  Research  Center 
IBM 

P.0.  Box  218 


84  07  24  014 

ISNJdXJ  lN3WNb3AO‘.>  IV  U3JlK10Hd3H 


REPRODUCED  at  government  fxpense 


Abstract 


An  explicit,  analytical  model  is  presented  of  finite  amplitude  waves  in 
shallow  water.  The  waves  in  question  have  two  independent  spatial  periods,  in 
two  independent  horizontal  directions.  Both  short-crested  and  long-crested 
waves  are  available  from  the  model.  Every  wave  pattern  is  an  exact  solution 
of  the  Kadomtsev-Petviashvili  equation,  and  is  based  on  a  Riemann  theta 
function  of  genus  2.  These  bi-periodic  waves  are  direct  generalizations  of 
the  well-known  (simply  periodic)  cnoidal  waves.  Just  as  cnoidal  waves  are 
often  used  as  one-dimensional  models  of  "typical"  nonlinear,  periodic  waves  in 
shallow  water,  these  bi-periodic  waves  may  be  considered  to  represent 
"typical"  nonlinear,  periodic  waves  in  shallow  water  without  the  assumption  of 
one-dimensionality. 
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1.  Introduction  , 

In  their  1895  paper,  Korteweg  and  deVries  showed  that  in  the  absence  of 
dissipation,  the  evolution  of  unidirectional,  long  water  waves  of  moderate 
amplitude  is  governed  approximately  by  an  equation  equivalent  to 


ft  +  6ffx  ♦  fxxx  -  0  • 


(1.1) 


This  has  become  the  standard  form  for  the  Korteweg-deVries  (KdV)  equation. 
They  also  discovered  periodic  solutions  of  (1.1)  in  the  form 


f ( x , t )  -  o2k2  en2  [o(x  -  ct  ♦  xQ);  k]  ♦  fQ  , 


(1.2) 


which  they  called  "cnoidal  waves".  In  (1.2),  cn[>;  k}  is  a  Jacobian  elliptic 
function  with  modulus  k  (0  <  k  <  1),  the  wave  speed  (c)  is  given  by 


c  -  6  iQ  -  1*  o2  (1  -  2k2)  , 


(1.3a) 


and  the  wavelength  (L)  is  given  by 


oL  -  2  K(k)  , 


(1.3b) 


where  K(k)  is  the  complete  elliptic  integral  of  the  first  kind  (cf.  Byrd  & 
Friedman,  1971).  Because  (1.1)  is  invariant  under  a  Galilean  transformation 
(t  -  t,  x  -  x  -  6  a-r,  f (x,  t)  -  v( x » t )  ♦  a),  we  may  normalize  its  periodic 
solutions  by  requiring 


Then  in  (1.2), 

io  '  "2“2  [|(kT  ‘  1  *  1,2  ]  • 

where  E(k)  is  the  complete  elliptic  integral  of  the  second  kind. 

Two  views  of  cnoidal  waves  are  popular.  Mathematicians  tend  to  regard 
the  KdV  equation  as  a  prototype  of  an  infinite  dimensional  Hamiltonian  system 
that  is  completely  integrable.  (For  the  analysis  of  (1.1)  with  periodic 
boundary  conditions,  see  Zabusky  &  Kruskal  (1965),  Novikov  (1974), 
Dubrovin  &  Novikov  (1974),  Dubrovin  (1975),  Lax  (1975),  Its  &  Matveev  (1975), 
Flaschka  (1975),  McKean  &  van  Moerbeke  (1975),  McKean  &  Trubowitz  (1976).)  Of 
particular  Importance  is  an  infinite  dimensional  family  of  exact, 
spatially-periodic  solutions  of  (1.1),  all  of  which  have  the  form 

f(x,t)  ■  232  in  0($i , . . . , ,  (1.6) 

where  ♦j  -  UjX  ♦  «jt  ♦  $j0,  and  6  is  a  Riemann  theta  function  of  genus  N.  One 
may  think  of  these  "genus  N"  solutions  of  (1.1)  as  the  periodic  analogues  of 
its  N-soliton  solutions,  which  can  be  obtained  from  (1.6)  in  an  appropriate 
limit.  In  the  simplest  case,  where  N  •  1,  the  solution  in  (1.6)  is  Identical 
with  that  in  (1.2). 

Ocean  engineers  often  adopt  a  different  viewpoint  that  was  popularized  by 
Wiegel  (I960).  From  this  perspective,  cnoidal  waves  are  important  not  as 
particular  solutions  of  a  partial  differential  equation,  but  rather  as 
one-dimensional,  analytical  models  of  "typical"  nonlinear,  periodic  waves  in 
shallow  water.  (In  this  paper,  we  use  "shallow-water  wave"  in  the  following 
sense.  Given  a  gravity  wave  with  amplitude  and  wavelength  L,  propagating 
in  water  of  depth  h,  its  Ursell  number  is 

U  -  L2/h3  .  (1.7a) 
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We  call  the  wave  a  shallow-water  wave  If 

U  -  0(1)  ;  (1.7b) 

see  Sarpkaya  &Isaacson  (1981,  p.215)  for  a  detailed  discussion  of  (1.7b).  The 
well-known  Stokes  expansion  of  water  waves  requires  U  <<  1 ,  and  becomes 
invalid  in  the  shallow  water  regime  of  interest  here.) 

The  practical  value  of  this  second  perspective  is  undeniable.  Jacobian 
elliptic  functions  are  well-documented  analytic  functions,  so  it  is  fairly 
easy  to  calculate  "typical"  values  of  wave  speeds,  forces,  mass  transport, 
etc.  for  engineering  purposes.  Cnoidal  wave  theory  assumes  that  wave 

amplitudes  are  not  too  large,  but  the  theory  has  been  generalized  by 

Laitone  (1961)  and  others  to  include  higher  order  nonlinear  effects.  Sarpkaya 
&  Isaacson  (1981)  summarize  these  generalizations,  and  give  a  concise 
description  of  cnoidal  wave  theory  as  a  practical  engineering  tool. 

One  of  the  major  shortcomings  of  cnoidal  wave  theory  as  a  practical  model 
of  shallow-water  waves  is  that  the  theory  is  one-dimen3ional,  whereas  the 
water  surface  is  two-dimensional.  It  follows  that  cnoidal  waves  are 
necessarily  long-crested,  whereas  both  long-crested  and  short-crested  waves 
are  observed  in  shallow  water. 

The  objective  of  this  paper  is  to  describe  a  two-dimensional 

generalization  of  cnoidal  waves.  The  waves  to  be  described  are  periodic  in 

two  independent  directions  on  the  water  surface,  and  they  are  proposed  here  as 

models  of  two-dimensional,  periodic  waves  of  finite  amplitude  in  shallow 
water.  (There  is  always  a  semantic  confusion  about  counting  dimensions  in 
these  problems.  Cnoidal  waves  have  two-dimensional  velocity  fields  and 

one-dimensional  surface  patterns;  in  this  paper,  we  call  cnoidal  waves 

one-dimensional.  The  waves  of  interest  here  have  three-dimensional  velocity 
fields,  and  two-dimensional  surface  patterns;  we  call  them  two-dimensional.) 
To  give  the  reader  some  notion  of  what  this  model  actually  is  and  in  what 
sense  it  generalizes  the  cnoidal  wave,  we  now  list  some  of  its  main  features. 
This  list  should  be  interpreted  as  a  series  of  claims,  which  will  be  verified 
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in  the  body  of  the  paper. 

a)  Just  as  the  KdV  equation  describes  the  evolution  of  long  water  waves  of 
finite  amplitude  if  they  are  strictly  one-dimensional,  the  KP  equation, 

(ft  +  6ffx  ♦  txxx  )  ,  +  3  fyy  *  0  .  (1.8) 

describes  their  evolution  if  they  are  weakly  two-dimensional  (Kadomtsev  & 
Petviashvili,  1970).  The  evolution  described  by  (1.8)  is  weakly 
nonlinear,  weakly  dispersive,  weakly  two-dimensional,  and  all  three 
effects  are  permitted  to  enter  at  the  same  order. 

b)  Both  the  KdV  and  KP  equations  are  Galilean  invariant,  so  their  periodic 
solutions  may  be  normalized  by  imposing  a  condition  like  (1.4). 

c)  Like  the  KdV  equation,  the  KP  equation  admits  an  infinite-dimensional 
family  of  exact  solutions  in  the  form 

f(x,y,t)  -  22 1  in  6  ,  (1.9) 

where  0  is  a  Riemann  theta  function  of  genus  N.  If  N  -  1 ,  these 
solutions  are  simply  cnoidal  waves. 

d)  The  waves  of  interest  in  this  paper  appear  at  N  -  2.  Then  the  KP 
equation  admits  solutions  in  the  form  (1.9)  that  are  periodic  in  two 
Independent  horizontal  directions.  Examples  of  these  two-dimensional 
waveforms  are  shown  in  Figures  2  and  3*  In  this  paper,  we  will  refer  to 
these  waves  as  "bi-periodic,"  or  as  "KP  solutions  of  genus  2."  (To  avoid 
confusion  we  emphasize  that  our  bi-periodic  waves  have  two  real  periods. 
They  should  not  be  confused  with  ordinary  elliptic  functions,  which  have 
two  complex  periods  but  only  one  real  period.) 

e)  The  cnoidal  wave  in  (1.2)  is  stationary  if  one  travels  with  the  speed 
defined  by  (1.3).  The  KP  solutions  of  genus  2  that  are  of  Interest  here 
also  are  stationary  if  one  travels  with  the  correct  (constant)  velocity. 
In  both  cases,  there  is  only  one  wave  pattern  for  all  time.  For  cnoidal 
waves  the  surface  pattern  is  one-dimensional;  for  bi-periodic  waves  it 
is  two-dimensional. 

f)  The  cnoidal  waves  in  (1.2)  have  three  arbitrary  parameters  (k,o,x0),  of 
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which  two  (k,o)  have  dynamical  significance.  These  two  parameters  can  be 
determined  uniquely  from  two  appropriate  measurements  of  the  wave.  The 
bi-periodic  waves  have  eight  arbitrary  parameters,  of  which  six  have 
dynamical  significance.  We  will  show  how  these  six  parameters  can  be 
determined  directly  from  a  few  measurements  of  the  wave, 
g)  The  derivations  that  lead  to  both  the  KdV  and  KP  equations  assume  that 
the  waves  in  question  are  too  small  to  break.  Neither  equation  ever 
predicts  breaking  waves,  but  either  equation  may  predict  waves  of  large 
amplitude  that  may  break  according  to  some  other  criterion. 

The  organization  of  this  paper  is  as  follows.  In  §  2  is  a  brier 
derivation  of  the  KP  equation  from  the  equations  for  inviscid  water  waves.  KP 
solutions  of  genus  1  (cnoidal  waves)  are  reviewed  in  §3.  The  KP  solutions  of 
genus  2  are  the  main  concern  of  this  paper,  and  these  are  described  in  detail 
in  §4.  How  to  use  these  solutions  as  models  of  bl-perlodlc  waves  in  shallow 
water  is  described  in  S5. 


2.  Physical  meaning  of  the  KP  equation 


The  KP  equation  is  a  two-dimensional  generalization  of  the  KdV  equation, 
and  its  derivation  as  an  approximate  model  for  water  waves  follows  closely 
that  of  the  KdV  equation.  Because  the  derivation  has  been  given  elsewhere 
(e.g.,  Ablowitz  &  Segur,  1979),  we  merely  review  the  main  points  here. 

In  the  simplest  case,  we  seek  the  irrotational  motion  of  an  inviscid, 
incompressible,  homogeneous  fluid,  subject  to  a  constant  (vertical) 
gravitational  force,  g.  We  employ  Cartesian  coordinates  (x>n>0  in  which  C 
increases  upward.  The  fluid  rests  on  a  horizontal,  impermeable  bed  at  (  *  -h 
and  has  a  free  surface  at  £  -  C(x»n,T).  Because  the  motion  is  irrotational  by 
assumption,  there  is  a  velocity  potential  *  that  satisfies 

V2*  -  0,  -h  <  £  <  c(x.n.T)  .  (2.1) 

It  is  subject  to  boundary  conditions  on  the  bottom  (£  »  -h), 

3^4>  »  0  ,  (2.2a) 

and  along  the  free  surface  (£  -  c(xtn.t)), 

iJ!{tMxVVn"*(  •  (2’2b) 

+  gC  +  —■  |V$|2  *  0  •  (2.2c) 

In  this  discussion,  we  are  interested  in  spatially  periodic  solutions  of 
(2.1),  (2.2),  and  we  require 

£im  /  dxc  (  X  +  X0.  n.T  )  -  0  ,  (2.3) 

This  is  a  normalization  of  c,  corresponding  to  (l.iJ).  It  implies  that  h  is  the 
mean  fluid  depth,  and  that  ;  »  0  in  the  absence  of  any  motion. 


Let  k  -  (fc,m)  be  a  representative  horizontal  wave  number  for  the  waves  in 
question,  with  **  -  t*  ♦  m*.  Orient  the  horizontal  coordinates  so  that  the 
X-direction  is  the  principal  direction  of  wave  propagation.  To  derive  the  KP 
equations  from  (£.1)  and  (2.2),  we  make  four  assumptions. 

(A)  Wave  amplitudes  are  small: 


(B)  The  water  is  rshallow  relative  to  typical  horizontal  wavelengths: 

(tch)2  «  1 

(C)  The  waves  are  nearly  one-dimensional: 

(m/£)2  «  1 

(D)  These  three  small  effects  all  have  comparable  influence: 

(m/£)2  -  0[(<h)23  -  0(e) 

When  these  assumptions  are  used  to  re-scale  the  original  equations 
systematically,  the  result  at  leading  order  is  a  one-dimensional  wave 
equation.  It  follows  that  the  surface  displacement  is  given  by 

C(Xtn,x)  -  eh  [  f  (x  -  *'ih  T;n,t)  *  f(x  +  /gh  r;n,T)]  ♦  0(e2h) 

where  f  and  F  are  known  in  terms  of  initial  data,  and  each  varies  slowly  in 
its  last  two  variables.  The  expansion  procedure  can  be  implemented  in  terms  of 


the  following  scaled,  dimensionless  variables, 


x  -  e1/2(x  ~  (gh) 1 /2i ]/h,  x  -  e1/2(x  ♦  (gh)1/2t)/h, 


(2.4 a) 
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y  -  en/h, 


t  -  e3/2(gh)1/2i/6h 


Then  we  look  for  solutions  of  (2.1),  (2.2)  of  the  form 


-  ~^[f(x,y,t)  ♦  F(x,y,t)j 


(2.^b) 


etc.  Secular  terms  arise  at  second-order  in  this  expansion  unless  the 
"right-going"  waves  satisfy  the  KP  equation, 

(ft  +  6f  fx  +  ^xxx^x  +  3fyy  “  ®  *  (1*8) 

and  the  "left -going"  waves  satisfy  an  equivalent  equation.  To  make  the  model 
quite  explicit,  we  also  write  the  KP  equation  for  the  right-going  waves  in  its 
dimensional  form  for  ?(x«n, t): 

[  <Bh)-'/2  t,  *  tx  *  f  i  c,  ♦  J-  Imjx  *  5  -  0  <2-5) 

The  components  of  fluid  velocity  in  the  (x»h,S)  directions  are,  for 
•h  S  5  S  C(Xth.t). 

u  -  (g/h)1/2  c  ,  (2.6a) 

v  -  (g/h )1/2/X  dx(3nC)  .  (2.6b) 

w  -  -(g/h)1/2(h  +  t)3x;  .  (2.6c) 

Note  that  only  the  vertical  component  of  velocity  (w)  has  any  vertical 
structure  at  leading  order.  To  the  same  order  of  approximation,  the  fluid 
pressure  below  the  free  surface  is  simply  hydrostatic: 

P  ~  pg(C  "  O  •  (2.6d) 
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Some  comments  about  this  derivation  of  the  KP  equation: 

(i)  Surface  tension  was  omitted  for  simplicity.  The  effect  of  including  it 
is  to  change  the  numerical  coefficient  of  Uxxxx)  in  (2.5).  For  water 
beneath  air,  this  change  is  negligible  if  the  water  depth  exceeds  a  few 
centimeters.  If  the  water  depth  is  less  than  about  1/2  cm,  however, 
this  coefficient  changes  sign  and  this  changes  the  qualitative 
character  of  the  solution  of  the  KP  equation.  We  recommend  caution 
when  using  the  model  described  here  for  periodic  waves  on  thin  films  of 
liquid. 

(ii)  This  derivation  neglects  dissipation,  which  is  assumed  to  be  even 
weaker  than  the  weak  nonlinearity,  weak  dispersion  and  weak 
two-dimensionality  included  in  the  KP  equation. 

(ill)  The  assumption  that  the  motion  is  irrotational  is  convenient  but  no* 
necessary,  as  shown  by  Benney  (1966). 

(iv)  The  original  equations  (2.1,  2.2)  are  invariant  under  a  horizontal 
rotation  of  coordinates.  This  symmetry  is  broken  in  choosing  the 
X-direction  as  the  principal  direction  of  wave  propagation.  Even  so, 
the  KP  equation  retains  a  remnant  of  this  symmetry;  (1.8)  is  invariant 
under  transformations  of  the  form: 

x  ♦  x  ♦  oy  -  3o2t,  3X  •»  3X, 

y  ♦  y  -  6ot,  3y  ♦  3y  ♦  ®3X,  (2.7) 

p 

t  *  t,  3^  •*  3t  -  6e3y  -  3«*  3X 

We  will  refer  to  these  transformations  loosely  as  "rotations,"  but  we 
emphasize  that  this  symmetry  group  only  approximates  actual  rotations 
in  the  physical  (x.n)  plane,  and  only  if  |3XJ  »  | | . 

Thus,  the  KP  equation  admits  waves  that  travel  in  arbitrary  directions 
in  the  (x,y)  plane,  but  we  can  expect  them  to  model  water  waves 


accurately  only  if  they  propagate  primarily  in  the  x-direction. 

(v)  The  KP  equation  also  exhibits  a  scaling  symmetry. 

x  ♦  Bx,  t  ♦  83t,  y  62y,  f  -*■  e_2f  .  (2.8) 

This  symmetry  is  inherited  from  the  arbitrariness  of  the  small 
parameter  (e)  in  (2.i»).  The  symmetry  parameters,  (a,$),  are  two  of  the 
six  parameters  available  in  the  model  of  doubly  periodic  waves  in 
shallow  water. 

(vi)  We  mentioned  below  (2.3)  that  (2.3)  defines  h  uniquely.  Because  of 

(2.6a),  (2.3)  also  Identifies  the  laboratory  coordinate  system  (i.e., 

X»n,t,T).  Integrating  (2.3)  over  /  dfc  shows  that  the  mean  horizontal 

J  -h 

motion  vanishes  in  the  x~directlon,  so  there  is  no  net  mass  transport 
in  that  direction  in  the  laboratory  frame. 


3»  KP  solutions  of  genus  1  -  cnoidal  waves 

Cnoidal  waves  have  been  discussed  at  length  elsewhere.  (For  a 
hydrodynamic  emphasis,  see  Wiegel,  I960  or  Sarpkaya  &  Isaccson,  1981;  for  a 
mathematical  emphasis,  see  Dubrovin,  1981  or  Boyd,  1982.)  We  review  them  here 
because  the  bi-periodic  waves  are  direct  generalizations  of  cnoidal  waves. 

3.1  Construction  of  KP  solutions  of  genus  1 

Real-valued  KP  solutions  of  genus  1  are  constructed  as  follows.  The 
phase  variable  is 

4>  -  yx  +  vy  ♦  wt  ♦  $0  ,  (3- la) 

or 

♦  -  u  (x  ♦  py  -  ct)  ;  (3.1b) 

♦  is  required  to  be  real.  A  theta  function  of  genus  1  is  defined  by  a  Fourier 
series: 

6  ( 4> I b )  -  exp('J  b  n2  ♦  in  4>j  (3-2) 

n— »  '  ’ 

0  is  both  convergent  and  real-valued  for  all  real  $  if  b  is  a  negative  real 
number,  which  we  now  require.  (Differences  in  notation  are  always  a  nuisance. 
0  ($)  and  b  are  called  t?  g  ($/2)  and  (2irit),  respectively,  by  Whittaker  and 
Watson,  1927.)  Krichever  (1976)  showed  that  the  KP  equation  (1.8)  has 


solutions  of  the  form 


f  (x,y,t)  -  23x2  In  6  (♦;b) 


(3.3) 


provided  c  is  a  particular  function  of  (y,  p,  b).  Standard  identities  for 
elliptic  functions  (cf.  Byrd  and  Friedman,  1971)  show  that  (3-3)  is 
equivalent  to 


where 


E(k) 

K(k)  ’ 


(3-4a) 


(3-^b) 


(3.4c) 


If  p  -  0  in  (3-D.  then  (3*4)  is  identical  with  (1.2)  and  (1.5),  the 
formulae  for  a  cnoidal  wave  solution  of  the  KdV  equation.  Its  speed  is  given 
by  (1.3a)  with  o  -  yK(k)/v.  For  p  *  0,  we  may  use  (2.7)  with  a  -  p  to  rotate 
this  KdV  solution  into  the  general  real-valued  KP  solution  of  genus  1,  with 
the  following  properties. 

i)  There  are  four  arbitrary  real  parameters  in  the  solution:  (y,  p,  $0)  in 
(3.1),  and  b  (b  <  0)  in  (3.2).  Of  these,  y  represents  the  scaling 
symmetry  in  (2.8),  p  represents  the  rotational  symmetry  in  (2.7),  and 
represents  the  translational  symmetry  of  any  of  the  independent 
variables  in  (1.8).  There  is  no  dynamical  significance  to  $0. 

11)  Given  these  parameters,  k  is  uniquely  determined  by  (3.4c),  the  waveform 
is  given  by  (3-3)  or  (3.4a),  and  the  speed  parameter  is  given  by 


(3.5) 


An  equivalent  expression  for  the  speed  in  terms  of  (y,  p,  b)  was  given 
by  Dubrovin  (1981). 

iii)  The  wavelength,  direction  of  propagation  and  total  speed  all  follow  from 
(3-D  and  elementary  geometric  considerations.  A  typical  cnoidal  wave 
is  shown  in  Figure  1. 


Instead  of  (3.2),  we  may  also  write 


©  (♦« b)  -(f£)1/2  X)  exp[t  (I " mit)2]  *  (3'6) 

This  representation  follows  from  an  identity  known  as  the  "Poisson  sum 
formula," 


^ ^ exp  (n2s  ♦  2niz  )  m  ~)1/2  ^  %xp  £  (  z  -  mir  )  y/sj 


(3.7) 


which  is  valid  for  Re(s)  <  0  (Whittaker  &  Watson,  1927,  p.  ^76).  For  any 
b  <  0,  it  follows  f row  (3.6)  that  the  two  terms  corresponding  to  m  •  0  and 
m  -  +1  are  larger  than  any  others  in  the  neighborhood  of  $  »  n.  If  b  *  0,  then 
in  this  region  these  two  terms  are  even  dominant  in  an  asymptotic  sense. 
[Similarly,  the  terms  corresponding  to  M  and  (M  ♦  1 )  dominate  in  the 
neighborhood  of  $  -  (2M+1)ir.  Thus  each  pair  of  terms  in  (3-6)  dominates  over 
exactly  one  period  of  6  (<J>),]  Retaining  only  the  terms  corresponding  to  m  -  0, 
1  in  (3.6)  and  substituting  into  (3-3)  yields  the  expression  for  single  KP 
soliton  riding  on  a  depressed  mean  level  (at  2 p  /b): 

f(x,y,t)~  222  lnfl  *  exp  -  x)jj*  2v2/b 

'  2(ff  Se0l,2[i(*  *  'j^TT  '  <3'e> 

This  is  correct  asymptotically  if  b  ■*  o  with  (p/b)  finite.  It  follows  that  the 
cnoidal  wave  can  be  viewed  as  a  superposition  of  overlapping  solitary  waves, 
placed  one  period  apart.  This  view  has  been  discussed  at  length  by  Boyd 
(1982). 


3.2  Reconstructing  a  cnoidal  wave  from  simple  measurements 

In  order  to  use  the  KP  solutions  of  genus  1  as  models  of  periodic  waves 
in  shallow  water,  one  must  be  able  to  infer  the  parameters  (u,  p,  b)  from 
measurements  of  the  wave.  This  can  be  done  in  a  variety  of  ways,  one  of  which 
is  outlined  here. 

The  cnoidal  waves  have  a  simple  characterization. 

i)  A  cnoidal  wave  is  real-valued  and  bounded,  for  all  real  (x,  y,  t). 
ii)  It  is  one-dimensional;  i.e.,  there  is  a  single  phase  variable 
£$  in  (3.1)],  which  is  linear  in  (x,  y,  t).  There  is  no  variation 
along  the  wave  crests  (i.e.,  along  $  -  const).  In  the  terminology  of 
oceanography,  cnoidal  waves  are  long-crested, 
ill)  A  cnoidal  wave  is  periodic  (in  $) . 

iv)  There  is  a  uniformly  translating  coordinate  system  in  which  the  cnoidal 
wave  is  stationary. 

Cnoidal  waves  are  the  only  KP  solutions  with  these  four  properties. 

Given  a  shallow-water  wave  with  these  four  properties,  one  can  determine 
the  parameters  of  the  corresponding  cnoidal  wave  in  three  steps,  as  follows. 

a)  Depending  on  the  problem,  one  may  or  may  not  be  able  to  orient  the 
X  -  axis  in  the  direction  of  wave  propagation.  If  so,  then  p  -  0  in 
(3*1)  and  in  (3.5).  If  not,  then  p  measures  the  angle  between  these 
two  directions;  one  can  expect  discrepancies  between  the  cnoidal  wave 
model  and  the  water  wave  if  p  is  large.  In  any  case,  p  can  be  measured 
directly  from  the  wave  pattern. 

b)  Once  p  is  known,  y  is  determined  from  p  and  L,  the  wavelength  in  the 
direction  of  wave  propagation: 

y  -  2tt[l2  (l  ♦  p2)  ]  "1/2  •  (3.9) 

c)  Given  y  and  p,  b  can  be  determined  from  the  wave  amplitude  or  from  the 
speed,  either  of  which  can  be  measured  directly.  Here  we  focus  on  the 
speed.  Suppose  y,  p  and  c  have  been  determined  by  direct  measurement. 
From  (3.5),  define 


(3.10) 


¥  •  -  -  "  jp2  •  (  — [ -3  —rr  ♦  2  -  k2 
p2  \  n  /  L  K(k) 

One  can  prove  that  y  is  a  monotonic  function  of  k,  and  also  of  b. 
Therefore,  b  is  determined  from  ¥.  Simple  approximations  are  obtained 
by  expanding  (3-4c)  and  (3*10)  for  k  *  1  and  for  k  ■»  0.  The  results  are 
as  follows.  As  b  *  -  •  (k  ■*  0), 


¥  -  -  1  ♦  24eb  . 


As  b  •*  0  (k  ■*  1), 


y 


For  -  •  <  b  <  0, 


-1 


♦  24eb  <  y 


+  3b  \ 

si' 


(3.na) 


(3.11b) 


(3.11c) 


Moreover,  one  may  approximate  ¥  by  (3.11a)  for  b  £  -  4.6.  and  by  (3.10b) 
for  b  ^  -  4.6,  with  an  error  that  never  exceeds  IS.  This  is  adequate 
for  comparisons  with  most  experimental  data. 
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4.  KP  solutions  of  genus  2 


The  KP  solutions  of  primary  interest  in  this  paper  have  genus  2.  They  are 
two-dimensional  generalizations  of  cnoidal  waves.  Aspects  of  these  solutions 
have  been  considered  by  Krichever  (1976),  Nakamura  (1979),  Krichever  and 
Novikov  (1980),  and  Dubrovin  (1981).  Related  work  by  Hirota  and  Ito  (1981), 
Bryant  (1982),  and  Finkel  and  Segur  (1984)  is  also  of  interest.  The  work  of 
Dubrovin  (1981)  is  particularly  important  from  our  perspective. 


4.1  Construction  of  theta  functions  of  genus  2 

The  construction  of  real-valued, nonvanishing  theta  functions  of  genus  2 

requires  a  certain  amount  of  machinery,  which  we  now  introduce.  The  first 

+  T 

ingredient  is  a  two-component  vector  of  phase  variables,  ♦  -  (fj,  $2)  ,  where 
♦j  “  4jX  ♦  vjy  ♦  «jt  ♦  #J0  ,  J  -  1,2  ,  (4.1a) 


-  Uj  (*  ♦  PjY  -  Cjt)  ♦  ♦jo  .  J  •  1,2  .  (4.1b) 

Both  components  are  required  to  be  real.  To  change  these  variables  into  those 
of  Dubrovin,  let  Wj  •*  iUj,  vj  ■*  iVj,  <dj  ■*  -41 W j ,  J  ■  1,2  . 


The  second  ingredient  is  a  real-valued  Riemann  matrix:  a  2x2,  symmetric 
matrix  that  is  real- valued  and  negative  definite.  Let 


be  real-valued,  g  is  negative  definite  if 

bn  <0,  b22  <  0,  bn  b22  -  b^2  >  0  .  (4.2b) 


1  £. 


A  theta  function  of  genus  2  is  defined  by  a  double  Fourier  series: 


m  m 

(♦l  •  *2’  fi  )  ■  ^  exp  j  j  m  •  g  *  o  ♦  i 


(4.3) 


where 


m  •  g  •  m  -  m|  b^  ♦  2  m^  m2  b^  ♦  bj 


m  *  $  ■  ♦  m2  ^2 

Clearly  e  is  real-valued.  It  converges  uniformly  in  $1 ,  <f>2  because  of 
( M . 2b ) .  It  is  periodic  in  ^  when  $2  *a  held  fixed,  and  vice-versa  (i.e.,  it 
is  quasi-per iodic) .  We  show  In  §  4.3  that  0  does  not  vanish. 

This  description  needs  further  refinement  for  two  reasons.  First,  two 
different  Riemann  matrices  can  be  equivalent,  in  the  sense  that  they  generate 
the  same  theta  function.  For  example,  let 


b  bX 


bX  bXf 


).*•(■  ). 

♦  d/  \  b( X+l)  t>(M)  ♦  d/ 


If  b  and  d  are  both  negative  and  X  is  real,  then  gj  and  g2  are  both  Riemann 
matrices*  Let 


in  m  ( f  in  2  J  9  n  -  [ ♦  m2),  m2~\ 


One  shows  by  direct  computation  that 


•  g2  *  01  ■  n  *  Bi  *  n 
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Then  it  follows  (by  renumbering  in  (4.3))  that  and  § 2  Produce  the  same 
theta  function,  and  are  equivalent.  Any  two  equivalent  Riemann  matrices  can 
be  related  to  each  other  by  a  symplectic  transformation  (see  Dubrovin,  1981, 
for  details). 

To  eliminate  this  ambiguity,  Finkel  and  Segur  (1984)  introduce  the  idea 
of  a  basic  Riemann  matrix. 

Definition:  A  real-valued,  2x2  Riemann  matrix  is  said  to  be  basic  if 


2  <  \j  *  2  ’  ^  “  1  ’ 2  * 


(4.4) 


where  Aj  -  bi2/bjj  •  They  also  prove  the  following  results: 

i)  Every  real-valued,  2x2  Riemann  matrix  is  equivalent  to  a  basic  Riemann 
matrix. 

ii)  Two  basic  Riemann  matrices  are  equivalent  only  if  their  A's  coincide. 

It  follows  that  every  real  theta  function  of  genus  2  can  be  generated  from  a 
basic  Riemann  matrix.  Henceforth,  we  consider  only  basic  Riemann  matrices, 
since  any  other  real  Riemann  matrix  can  be  reduced  to  its  basic  form  via  an 
algorithm  given  by  Finkel  &  Segur  (1984). 

Basic  Riemann  matrices  have  a  natural  representation  as  follows.  Define 


b  -  max  ^>1 1 ,  b22^ 


A  -  b12/b 


d  -  det  g  /b 


(4.5) 


Then  every  basic  Riemann  matrix  takes  one  of  two  forms 


(b  bX  \  /bX2*d  bx\ 

)  or  g2  -[  )•  <**6a> 

bx  bX^+d  /  \bX  b  f 

With 

b<0,  l2Sj,  d$b(  1-X2  )  .  (4.6b) 

The  relation  between  these  two  forms  is 

©  (♦!»  ♦2*  6l)  “  ©  ^2*  ©i »  82 )  *  (4.6c) 

For  definiteness,  we  will  always  use  g1  unless  otherwise  noted. 

The  second  difficulty  with  (4.3)  is  that  it  admits  degenerate  cases.  A 
2x2  Riemann  matrix  is  called  decomposable  it  is  is  equivalent  to  a  diagonal 
matrix,  and  Indecomposable  otherwise.  For  example,  in  (4.2)  £  is  decomposable 
if  b12  -  0.  In  this  case  0  in  (4.3)  becomes  a  product  of  two  theta  functions 
of  genus  1.  Dubrovln  (1981)  shows  that  theta  functions  of  genus  2  generate 
nontrivial  KP  solutions  if  and  only  if  they  are  indecomposable.  He  also  gives 
an  explicit  test  for  decomposability,  which  requires  the  evaluation  of  a  4x4 
determinant.  A  simpler  test  is  given  by  the  following  Theorem  (Finkel  and 
Segur,  1984):  Let  g  be  a  basic  Riemann  matrix,  as  in  (4.6).  Then  g  is 
decomposable  if  and  only  if  X  ■  0. 

Summary : 

Let  (b,  d,  X)  be  three  real  parameters  satisfying  (4.6b),  with  X  4  0. 
These  parameters  generate  a  basic  Riemann  matrix  in  the  form  (4.6a).  Every 
real,  2x2,  indecomposable  Riemann  matrix  is  equivalent  to  a  basic  Riemann 
matrix.  Basic  Riemann  matrices  with  X  *  0  generate  real-valued,  indecomposable 
theta  functions.  These  can  be  written  in  the  form 


4.2  Construction  of  KP  solutions  of  genus  2 

Given  an  indecomposable  theta  function  of  genus  2,  KP  solutions  are 
constructed  in  the  usual  way: 

f(x,y.t)  -  2  £n  ©  (  8  )  • 

Dubrovin  (1981)  proves  that  f(x,y,t)  actually  satisfies  the  KP  equation, 
provided  that  its  parameters  satisfy  certain  relations.  To  express  these 
relations,  we  need  two  more  concepts:  theta-constants,  and  two  more  phase 
variables. 

Let  p  be  a  two  component  vector,  which  can  take  on  one  of  four  values: 

»•(;)•(:)■(’? )•(-)•(“«)• 

Every  Riemann  matrix  generates  a  four-component  theta  constant  of  the  form 


5  tn-  £  £ 

m,  m2  -  — 


(4.10a) 


where  m  •  (m^,m2),  and  the  product  is  defined  below  (4.3).  If  8  Is  *n  basic 
form  (for  definiteness,  in  (4.6a))  then  this  can  be  written  as 


©  [  P  ]  *^  exp  d  (  m2  ♦  p2  )  2  £  exp  j  b  (  m,  ♦  p,  ♦  A  (m2  ♦  P 2)f  • 

m2  ’  '  m1  I  (4.10b) 

-  .* 

Obviously,  6  [p]  is  differentiable  with  respect  to  (b,  d  and  A),  and  each  of 
these  derivatives  is  a  4-component  vector  Indexed  by  p.  Define 

♦l|  -  #2  “  *♦]  .  $3  -  ^i  -  A  4.2  »  (4.11a) 
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where 


i 


i  - 


bx 


J1 2 


b22  bx2  ♦  d 


(it. lib) 


We  show  in  §  4.3  that  the  wave  crests  lie  on  +3  •  const  and  on  44  «  const.  The 
corresponding  wave  numbers  are 


H3  -  b. 


-  X 


w2 


bjj  "  »*2  “  *  **1 


(4.11c) 


with  similar  definitions  for  (v^,  vjj,  u^). 

Now  we  may  state  the  main  Theorem  (Dubrovin,  1981):  Let  g  be  a  real. 
Indecomposable  Riemann  matrix;  i.e.,  satisfying  (4.6)  with  X  4  0.  Then  a 
function  of  the  form  (4.8)  solves  the  KP  equation  (1.8)  if  and  only  if  its 
parameters  are  related  by  the  four  equations.  Indexed  by  p,  in  (4.12). 


M  t  -  4  §  $ 


where 


t  -  I  bj  u1  +  3v2  ^ 

b^  ujj  4  b|j  w1  6v1  vjj 


b|j  ♦  3v^ 


V 


“I 


4b|  bjj 


2  2 
6b‘  b$ 


4h  bj 


b{ 


^  ■  3^  0[p],  ~  3^  0[p],  3(3  0[p],  0[p]  j 
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(4.12) 


'  %' 


~  1 


1 

•  J 

*  ’-\M 


and  D  is  a  constant  of  integration  with  no  physical  significance.  H  is 
invertible  if  and  only  if  g  is  indecomposable,  so  (4.12)  always  can  be  solved 
for  t. 


After  is  inverted  in  (4.12),  the  first  and  third  of  the  resulting 
equations  can  be  written  as 


(4.13a) 


2  2 
Cj,  -  3  Pj|  UJ4  <*>1,  *  3 


Here  P^(z;4)  and  Pj,(z;4)  are  well-defined  polynomials  of  degree  4.  Then 
eliminating  (D,  u1 ,  <uj,)  from  (4.12)  yields  one  more  relation. 


(4.130 


where  P(z;6)  is  a  well-defined  polynomial  of  degree  6.  The  resulting  solution 
will  be  real-valued  if  ujj/ui  is  chosen  so  that 


(4.13d) 


To  summarize,  (4.8)  provides  a  real-valued  solution  of  the  KP  equation, 
(1.8),  whenever:  (i)  g  is  indecomposable  (i.e.,  A  ^  0  in  (4.6));  and  (li) 
(4.13)  is  satisfied.  There  are  no  other  constraints. 
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The  conditions  in  (4.13)  can  be  interpreted  in  terms  of  the  wave  pattern 
as  follows: 


(i)  According  to  (4.11c), 


(4.14) 


(ill  If  pj  v2  -  u2  vi  •  then  at  any  fixed  t,  ^  -  constant  and  $2  *  constant 
on  parallel  lines  in  the  (x,y)-plane.  Using  (2.7),  such  a  solution  can  be 
"rotated"  so  that  -  v2  -  0;  i.e.,  these  are  actually  KdV  solutions  of 
genus  2.  According  to  (4.13c)  and  (4.14),  pjv2  *  u2v1  wherever 
P  0.  In  this  way,  every  indecomposable  2x2  Riemann  matrix 
generates  exactly  6  KdV  solutions  of  genus  2,  corresponding  to  the  6  roots  of 
P  ^ u^/ui ; 6 J  -  0.  The  solutions  may  or  may  not  be  real.  This  identification 
is  due  to  Dubrovln  (1981).  Our  interest  in  the  present  paper  is  in 
two-dimensional  waves,  so  we  will  systematically  neglect  these  one-dimensional 
KdV  solutions. 

(ill)  Suppose  u1v2  4  y2vi •  Then 


//  d^1  /\  d$2  -  (  u1v2  -  u2vl) J* f dx  ^  dy  (4.15) 

represents  a  non-zero  area-element  in  the  (x,y)-  plane.  For  any  fixed 
("♦l#  "?2  )*  a  period  parallelogram  of  9  (  ^ ,  $2;§)  has  as  its  vertices  four 
points  defined  by 


(♦v  *2  )  "  (  *1  •  *2  )  •  (  *1  +  2*»  *2  )  .  (  ♦  2n,  $2  ♦  2v  )  ,  (  ,  $2  ♦  2ir  ) 


The  area,  A,  of  this  parallelogram  is  given  by 


so  that 


(4.16) 


It  follows  that  (^ .  1 3c)  can  be  interpreted  as  specifying  the  area  of  the 
period  parallelogram  of  each  KP  solution.  As  one  might  guess,  this  area  is 
Independent  of  rotations  (2.7),  but  depends  on  scaling  (2.8). 


(iv)  If  one  seeks  a  "nonlinear  dispersion  relation"  for  these  waves  of 
genus  2,  it  is  given  by  (4.13a,b).  These  two  relations  are  the  generalization 
to  genus  2  of  (3.10).  and  ¥4  are  independent  of  either  of  the  symmetries  in 
(2.7)  or  (2.8).  They  depend  only  on  g  and  on  (  They  determine  the  wave 
speeds. 


(v)  Suppose  u-|V2  4  p2v1  ’  Then  at  any  $1  “  constant  and  $2  -  constant 
intersect  transversely  at  a  unique  point  (x(t),  y ( t ) ) .  As  t  changes,  this 
point  moves  in  the  (x,y)-plane  with  a  constant  speed  given  by 

—  .  J_£ - £_L  ,  ( 14 . 1 7a ) 

dt  -  p2v1 


d L  m  ~  ^1  u2 

dt  p']  v2  -  p2Vj 


(4.17b) 


where  the  parameters  are  determined  by  (4.13).  In  fact,  every  point  in  the 
wave  pattern  moves  with  this  same  speed.  The  entire  wave  pattern  is  stationary 
in  a  coordinate  system  that  translates  uniformly  with  a  velocity  given  by 
(4.17).  Every  KP  solution  of  genus  2  with  y1v2  *  p2v1  is  stationary  in  a 
coordinate  system  appropriate  to  that  wave,  just  -as  every  cnoidal  wave  is 
stationary  in  an  appropriately  translating,  one-dimensional  coordinate  system. 


For  completeness,  we  now  summarize  a  different  description  of  these  KP 
solutions  of  genus  2,  in  terms  of  the  underlying  algebraic  geometry.  More 
details  are  given  by  Dubrovin  (1981).  Every  indecomposable  2x2  Riemann  matrix 
(real  or  not)  corresponds  to  a  compact  Riemann  surface  of  genus  2.  This 
complex  surface  is  a  hyperelliptic  curve,  defined  by  an  equation  of  the  form 

w2  -  P  (Z;  5)  or  w2  -  P  (Z;  6),  (4.18) 

where  P  (Z;  n)  is  a  polynomial  of  degree  n.  The  surface  is  topologically 
equivalent  to  a  sphere  with  two  handles.  Given  such  a  surface,  there  is  a 
well-defined  procedure  to  generate  a  Riemann  matrix  on  the  Riemann  surface. 
Moreover,  given  a  surface  and  a  fixed  reference  point  on  it,  there  is  a 
mapping  (the  Abel  mapping)  from  each  point  on  the  Riemann  surface  to  a  point 
on  a  complex  2-torus,  called  the  Jacobi  variety.  This  torus  is  coordinatized 
by  two  independent  complex  variables.  For  example,  if  we  allowed  and  $2  in 
(4.1)  to  take  on  complex  values,  they  would  coordinatize  the  torus.  Each  of 
these  complex  variables  has  two  complex  periods,  and  at  most  one  real  period. 

The  inverse  problem  may  be  stated  as  follows.  Given  a  Riemann  matrix,  g, 
did  it  come  from  a  Riemann  surface?  If  so,  from  which  one?  For  genus  2,  the 
solution  as  follows:  g  corresponds  to  a  Riemann  surface  of  genus  2  if  and 
only  if  g  is  indecomposable.  Here  is  how  to  construct  the  curve  (4.18) 

explicitly,  if  g  is  indecomposable.  Construct  a  theta  function  as  in  (4.3), 

and  use  (2.8)  to  scale  u-j  -  1 .  Then  (4.13c)  has  the  form  (4.18),  and  defines 
the  hyperelliptic  curve.  Points  on  the  Riemann  surface  are  identified  by  the 
complex  variable  (u^/u^)'  This  surface  has  6  branch  points  (called 
Welerstrass  points),  defined  by  P  (Z;  6)  -  0.  These  6  points  represent  the 
6  KdV  solutions  discussed  above. 

Every  KP  solution  of  genus  2  has  eight  free  parameters:  (b,d,A)  in 

(4.6),  and  ( ,  u2 ,  ,  4^  q  ,  $2q  ^  ^  (4.1).  Then  ( v2 ,  ,  u2 )  in  (4.1)  are 

determined  by  (4.13).  Of  the  eight  parameters,  and  $20  affect  neither  the 
dynamics  nor  the  algebraic  geometry.  Of  the  remaining  six,  we  may  always 


choose 


Ml  -  1  and  -  0  , 

by  using  (2.7)  and  (2.8).  Thus  each  KP  solution  of  genus  2  has  just  four 
parameters  that  cannot  be  changed  by  a  Lie  point-symmetry  of  (1.8).  Those  are 
(b,d,X),  which  determine  the  underlying  Riemann  surface,  and  which 
determines  the  fixed  reference  point  on  that  surface. 
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4.3  A  floliton  representation 


By  using  twice  the  Poisson  sum  formula  (3.7).  one  easily  shows  the 
equivalence  of  the  following  three  representations: 


(4.7) 


(4.19a) 


where  all  sums  are  taken  over  (-•,»).  It  is  numerically  efficient  to  use 
(4.7)  if  both  b  and  dare  large  (and  negative),  to  use  (4.19a)  if  d  is  large 
with  b  small,  and  to  use  (4.19b)  if  both  b  and  d  are  small.  (The  other 
possibility,  b  large  with  d  small,  is  excluded  by  (4.6).)  Note  from  (4.19b) 
that  if  (b,  d.  A,  are  a11  r®al,  then  e  is  the  sum  of  positive  terms. 

Therefore  0  never  vanishes,  and  the  associated  KP  solution  is  real  and 
bounded.  The  corresponding  formulae  for  the  4-component  vector  of  theta 
constants  are: 


©  [p]  -  ]C  exp  j  d  ( m2  ♦  P2  )  2  }  * 

£  exp  |  b  (  m,  ♦  p,  ♦  A(m2  P^))2}  » 

St53.(^)’/2  exp  j  a  (®s  ♦  p2)2}  •[ 


(4. 10) 


♦  2  V.  exp 
n=1 


|(-r)2| *cos  {2nni  (pi +  x  (m2 +  p2))|J* 


(4.20a) 


0  [p] 


(bd  ) 


j,2 


2  1 


2 


-  2iriri2P2 


(4.20b) 


We  now  propose  a  soliton  interpretation  of  the  KP  solutions  of  genus  2. 

a)  These  solutions  can  be  viewed  as  a  doubly  periodic  array  of 

interacting  solitons.  This  description  becomes  exact  in  the  limit  b-»o,  d-»o. 

b)  Away  from  the  regions  of  interaction,  the  crests  of  the  individual 

solitons  lie  along  $2“const  and  ^-const,  as  defined  in  (4.11).  (Within  a 

region  of  interaction,  of  course,  individual  solitons  cannot  even  be  defined 
unambiguously.)  Thus  the  edges  of  the  period  parallelogram  lie  along  (f1 -const, 
and  $2*const*»  but  t*ie  wave  crests  ordinarily  do  not.  According  to  (4.11), 
therefore,  X  and  A  determine  the  angles  of  rotation  of  the  wave  crests  from 
the  directions  of  periodicity.  This  is  a  direct  geometric-dynamical 

interpretation  of  the  off-diagonal  term  in  the  Riemann  matrix:  it  determines 
how  strongly  the  waves  of  the  two  families  ($2  and  $4)  affect  each  other. 


-V  '*  _*«  *•  ■  •  ,*•  .  •  .**  «  *  «*•  ‘  *  *  V*  V 


-  \ 


c)  A  given  wave  crest  experiences  a  phase  shift  from  each  of  its 
interactions  with  waves  of  the  other  family.  The  strength  of  the  phase  shift 
is  proportional  by  X  (or  X). 

The  soliton  interpretation  can  be  proved  in  the  limit  b+0,  d-*0,  so  we 
consider  that  case  first.  Note  that  (4.19b)  can  be  written  as 


(4.21) 


In  the  limit  of  b-K),  d+0  and  in  a  neighborhood  of  (^  -  »,$2  -  v),  the  four 
terms  in  (4.21)  corresponding  to  (nj-  -1,  0),  (n2-  “1,  0)  dominate.  Retaining 
only  those  terms  in  (4.21)  yields 


(bd) 


exp 


exp 


exp 


(4.22) 
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Substituting  this  into  (4.3)  yields  the  approximate  KP  solution 


f  (  x,y,t)  -  232  ln£l  ♦  exp  1^3  ♦  exp  ni|  ♦  exp  ^3+  +  AjJ 


♦ 


2y2  2  (u2  “  ^ )2 

~b~  *  d 


(4.23) 


where 


n4  *  "  “  ( *2  ‘  A  "  ¥)  *  “  ~  (  *4  "  11 )  * 

-  2lt(i  4  t)  (*1  '  ’)  4  2XT  *2/d  -  -  4  V)  (*3  “  w)  • 


4ir2A 


and  we  have  used  (4.6)  and  (4.11).  The  parameters  in  (4.23)  are  constrained  by 
(4.13),  which  in  this  limit  becomes 


(4.24) 


2iryJj^  ♦  d)  ^ 

2  -  0 

,  d  /  \  bd  /  . 

\p4  p3/ 
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2  _(p  _  p  \2 
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Note  that  the  wave  -speeds  remain  finite  in  this  limit  only  if  (uij/d)  remains 
finite  as  d-*0,  and  if  (i^/b)  remains  finite  as  b+0.  But  (4. 23)  and  (4.2*0 
represent  an  exact  solution  of  (1.8)!  They  describe  two  solitons,  interacting 
obliquely  on  a  depressed  mean  level  (cf.  Satsuma  (1976)  or  Ablowitz  &  Segur 
1981,  p.  189).  Away  from  the  region  of  Interaction,  the  crests  of  the 
individual  solitons  lie  along  $3  ■  constant  and  ■  constant.  The  phase 
shift  of  the  -  soliton  due  to  the  interaction  is  obtained  from  (4.23)  by 
letting  $3  ■*  ±  •,  holding  ^  fixed.  The  result  is 

A  ■  2vl  ,  A  43  m  2x1  •  (4.25) 

The  largest  neighborhood  of  (^  ■  x-j  $2  *  *)  in  which  the  four  terms  in 
(4.22)  dominate  is  precisely  a  period  parallelogram  centered  at 
(♦i  -  i.  -  *)•  Moreover,  in  the  corresponding  period  parallelogram 
centered  at  (+,  -  (2  N,  ♦  1)x,  $2  -  (2N2  +  1  )  there  are  ln  t1**21)  exactly 
four  dominant  terms,  indexed  by  (n1  ♦  Nj  -  -  1,  0),  (n2  ♦  N2  ■  “1,  0).  Thus 

the  entire  series  in  (4.21)  represents  a  bi-periodic  repetition  of  a  2-soliton 
solution  of  (1 .8). 

To  summarize,  we  have  now  established  the  validity  of  the 
soliton-interpretation  given  above  in  the  limit  b  0,  d  ♦  0.  Figure  2a  shows 
a  two-soli ton  solution  of  (1.8).  Figure  2b  shows  a  genus  2  solution  with 
(b  -  -  1 ,  bl2  ♦  d  -  -  1).  The  sense  in  which  the  solution  of  genus  2 
represents  a  bi-periodic  extension  of  the  two  soliton  solution  is  evident. 

What  happens  if  b  and  d  are  not  small?  Given  any  b,d  satisfying  (4.6), 
then  in  a  period  parallelogram  centered  at  ($j  ■  x,  $2  -  x),  the  four  terms  in 
(4.22)  are  bigger  than  any  of  the  others  in  (4.21),  but  they  are  not 
necessarily  dominant  in  an  asymptotic  sense.  Even  so,  there  is  a  practical 
sense  in  which  the  wave  crests  still  lie  along  $3  -  constant  and 
♦jj  ■  constant,  provided  this  statement  is  interpreted  correctly.  This 
assertion  is  corroborated  in  Figure  3,  which  shows  some  typical  KP  solutions 
of  genus  2  outside  the  soliton  limit.  Each  wave  pattern  consists  of  a 
bi-periodic  array  of  peaks,  connected  to  each  other  in  two  directions  by 


ridges.  The  peaks  are  the  interaction  regions,  where  the  two  wave  amplitudes 
"add".  If  those  peaks  were  well  separated  as  in  Figure  2,  the  ridges  would 
become  the  wave  crests  of  the  individual  solitons.  In  Figure  3  the  peaks  are 
not  well  separated,  and  the  "individual  wave  crests"  are  simply  the  lines  of 
steepest  descent  and  ascent  across  the  saddles.  We  have  sketched  lines  of 
constant  43  and  constant  +4  in  Figures  2c, 3d, and  3g  to  show  that  they  identify 
the  wave  crests  (i.e.,  the  lines  of  steepest  ascent  in  the  saddles)  in  the 
general  case. 

Do  the  wave  crests  still  lie  along  $3  -  constant  and  4^  -  constant,  even 
if  b,d  ♦  -  •? 

(i)  Suppose  d  *  -  •,  with  b,  A  finite.  Then  A  ♦  0  from  (4.11),  so 
♦3  ■*  4^.  In  this  limit,  therefore,  the  assertion  is  that  the  crests  lie  along 
4j  -  constant  and  4j,  -  constant.  Its  validity  for  d  -  •  follows  easily  from 
(4.19a);  it  is  shown  in  Figure  1,  which  is  actually  a  wave  of  genus  2  with 
b  -  -  3,  d  -  -  12,  A  -  0.3,  -  1*2  “  °*5, 

(ii)  Suppose  b  -»  -  ■  and  d  ■*  -  •.  Then  the  appropriate  representation  of 

0  is  (4.7).  We  will  show  in  §4.4  that  the  limit  b  ♦  -  •,  d  ♦  -  •  is 

«■» 

nondegenerate  only  if  A  *  0  as  well.  In  this  limit,  therefore,  A  ♦  0,  A  ■*  0 
and  the  wave  crests  lie  along  the  directions  of  periodicity. 

In  summary,  we  assert  that  the  wave  crests  always  lie  along  43  ■  constant 
and  4j|  •  constant.  We  have  proved  this  assertion  in  certain  limiting  cases, 
and  have  offered  supporting  evidence  in  other  cases. 

Now  we  turn  to  point  (c)  above.  We  have  already  shown  that  in  the  limit 

b  ♦  0,  d  ♦  0,  the  phase  shifts  in  each  period  of  the  lines  43  •  constant  and 

4jj  ■  constant  are  given  by  (4.25).  In  fact,  the  validity  of  (4.25)  does  not 
rely  on  any  special  limits;  it  is  a  consequence  of  elementary  geometry.  In 
Appendix  A,  we  prove  the  following  theorem.  Let  4;  and  42  be  independent  real 
variables.  For  fixed  A  and  A  (0  &  A^  S  A^  S  1/4),  define  43  and  4ij  by  (4.11) 
on  0  i  41  <  2 it,  0  i  42  <  2v,  and  by  periodic  extension  outside  this  rectangle. 
Then  once  each  period,  43  -  constant  and  4^  -  constant  experience 
discontinuous  shifts  whose  size  and  direction  are  given  by  (4.25). 


4.4  Limits  and  asymptotic  formulae 

The  KP  solutions  of  genus  2  generalize  and  unify  several  simple  models  of 
waves  in  shallow  waters,  including  one-dimensional  cnoidal  waves,  obliauelv 
superposed  infinitesimal  waves  (i.e.,  Fourier  modes),  and  2-soliton  formulae. 
He  show  next  how  these  previously  known  results  emerge  from  the  KP  solutions 
of  genus  2  in  various  limits.  He  also  develop  asymptotic  formulae  that  will  be 
useful  in  §5. 

(a)  Limit  of  one  weak  wave:  d  • 

If  d  •  in  (4.7),  then 


-  ]E  exp  1 1  b  a2  ♦  1  Mi  I 

■  ) 


♦  e  d/2  exp 


exp  ^  i  X*  exp  j  j  b  ^m  «•  X  j  ♦  i  (n  ♦  a)  ^  j  ♦  (*) 


(•“)  • 


(4.27) 


where  (*)  denotes  complex  conjugate.  Comparing  this  with  (3*2)  shows  that  at 
leading  order,  (4.27)  represents  a  one-dimensional  cnoidal  wave  with  phase 
variable  ^ .  There  is  no  phase  shift  of  the  cnoidal  wave  at  this  order  because 
the  second  wave  is  too  weak.  The  next  order  describes  a  weak  perturbation  of 
the  cnoidal  wave,  with  an  amplitude  o(ed/2)  .  This  perturbation  solves  the 
equation  obtained  by  linearizing  (1.8)  about  a  cnoidal  wave.  He  will  refer  to 
these  weak  perturbations  as  "linear  waves".  An  example  of  such  a  perturbed 
cnoidal  wave  is  shown  in  Figure  3a. 


The  limiting  form  of  (4.12)  depends  on  whether  or  not  b  remains  finite  as 
d  -*  Suppose  b  remains  finite.  Then  for  either  p1  •  0  or  p1  ■  1/2  in 

(4.10), 

eCo1  ]  -  52  exp  jb^m  +  pj)  2  j  ♦ 

ed  52  \  exp  |  b  (  m  ♦  p1  ♦  A  \  2  j  *  exp  |  b  ^  m  ♦  p1  -  A 


-  0  (  e*d  )  . 

6[?/i]  •  e  d/1*^52  exp  |  b  (  m  +  p1  ♦  A/2  J  2  j 

♦  exp  |  b  ( m  ♦  p1  -  A/2  j  2  j  ♦  0  ( e2d )  j . 

When  these  limiting  forms  are  substitute  into  (4.12),  the  two  equations 
with  p2  -  0  become 

D  6,  [o1]  -  (uj  w,  ♦  3  v,2)  3b  ©t  [o1]  ♦  4y{  0,  [o1]  -  0,  (4.28) 


where 


6,  [  ?’]  -  E 

m 


exp 


{  b  (m  +  p, 


These  are  exactly  the  equations  one  obtains  for  a  cnoidal  wave.  If  D  is 
eliminated  from  (4.28),  the  resulting  equation  is  equivalent  to  (3.10).  The 
other  two  equations  in  (4.12),  with  p2  ■  1/2,  determine  the  speed  and 
orientation  of  the  linear  perturbation. 


(b)  The  small  amplitude  limit:  d  +  -  •,  b  ■*  -  •,  yj  -  0  (1) 


In  this  limit,  (4.7)  becomes 


corresponding  to  an  approximate  KP  solution  of 


H 


f  (  x,y,t )  -  2pf  exp  |  j  |  cos  ^  +  2p2  exp  j  ~  (bA2  +  d  )  |  cos  $2  , 


(4.29) 


where  and  $2  are  defined  by  (4.1).  Clearly  (4.29)  represents  a  linear 
superposition  of  two  Fourier  modes  of  small  amplitude,  and  one  expects  (4.13) 
to  produce  the  dispersion  relation  of  the  linearized  problem  in  this  limit. 
This  expectation  is  correct.  One  shows  that  in  this  limit,  (4.13)  becomes 
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F 

(4.30c) 


where  we  have  used  (4.1b).  Because  b  ■*  -  »  with  ,  p2  finite,  (4.30c)  can  be 
satisfied  only  if  A  +  0  as  well.  Then  one  may  prescribe  two  Fourier  modes  in 
(4.29)  arbitrarily,  and  (4.30c)  determines  (bA).  Because  A  +  0,  the  wave 
crests  and  the  lines  of  periodicity  coalesce  in  this  limit. 

For  use  in  §5,  we  also  give  the  next  corrections  to  (4.30a,b).  As 
b  -  »,  d  +  -  ®,  with  (bA)  finite, 


"  3  Pj  /  p  \ 

-  -  — - —  -  -  1  +  24  exp  (b)  -  12  exp  (  bA*  ♦  d)  •  M/ 
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-1+24  exp  (  bA2  +  d)  -  12  exp  (b)-  M/p26  , 


(4.3D 
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v2  [  (“1  *  ^2)  2  exP(b^)  ♦  (  V]  ~  V2  )2  exp(-bA)  -2  ^  My2  *  w22  )  j 


Cc)  Limit  of  one  soliton  plus  perturbation:  d  b  0 

Consider  (4.19a)  in  the  limit  b  0,  d  ■»  -  ®,  and  in  a  neighborhood  of 


e '(^)'/2  exp  I  ’  * exp  It  (' •  ♦>  )l  *  l"-32) 

+  2  exp  ^  ^  cos  4>ij  ♦  exp  |  ^  (  "  ~  $1 )  j  cos  (  4>l|  +  2*x)^ 

where  $.j  and  4^  are  defined  in  (4.1),  (4.11).  The  parameters  of  this  solution 
satisfy 


(4.33a) 


2  J| 

+  3vjj  -  p^  -  0 


(4.33b) 
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(4.330 


At  leading  order  in  exp  (d/2)  ,  (4.32)  gives  a  one-dimensional  soliton  as 

b  *  0.  Equation  (4.33a)  shows  that  the  speed  of  this  soliton  remains  finite  as 


b  -*  0  only  if  ij]  +  0  so  that  (wj/b)  remains  finite.  At  second  order,  (4.33b) 
shows  that  the  speed  of  the  linear  wave  is  unaffected  by  the  presence  of  the 
soliton  of  leading  order.  The  next  order  corrections  to  (4.33a,b)  are 
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(d)  Limit  of  two  solitons:  b  -*•  0,  d  0 


This  limit  was  already  discussed  in  i  4 . 3 »  and  from  quite  a  different 
viewpoint  by  McKean  (1979).  The  higher  order  corrections  to  the  wave  speeds 
are 


(4.35) 


Note  that  both  ¥3  and  are  singular  in  this  limit,  although  the  wave  speeds 


5.  Using  wave  measurements  to  calibrate  the  model 

The  KP  solutions  of  genus  2  have  eight  free  parameters:  b,  A,  d  in 
(4.6),  and  ^ ,  y2,  v1 ,  $10,  #20  in  (4.1).  In  order  to  use  these  solutions  as 
models  of  two-dimensional  periodic  waves  in  shallow  water,  one  must  be  able  to 
infer  these  parameters  from  measurements  of  the  water  wave.  The  purpose  of 
this  section  is  to  show  how  to  infer  these  parameters  frco  two  different  kinds 
of  wave  measurements.  At  this  time,  we  have  not  proven  rigorously  the 
validity  of  every  step  of  our  algorithms,  but  we  have  identified  the  major 
unproven  points  in  five  conjectures. 

Certainly  the  KP  solutions  of  genus  2  cannot  describe  all  possible  waves 
in  shallow  water.  Every  two-dimensional  KP  solution  of  genus  2  that  is 
generated  by  a  real  Riemann  matrix,  as  discussed  in  i  4,  necessarily  has  the 
following  four  properties,  so  the  water  waves  that  they  model  should  also  have 
these  properties. 

(i)  The  solution  is  real-valued  and  bounded,  for  all  real  (x,  y,  t). 

<ii)  It  is  spatially  periodic  in  two  independent  directions,  which  we  may 
identify  with  ^  and  $2  in  (4.1).  If  ^  and  $2  are  not  collinear,  then 
W1V2  *  W2V1  ’  We  call  a  solution  two-dimensional  if  *  y2vi  • 

(iii)  There  is  a  uniformly  translating  coordinate  system  in  which  the  entire 
two-dimensional  wave  pattern  is  stationary. 

(iv)  The  common  features  of  two-dimensional  KP  solutions  of  genus  2  can  be 
inferred  from  the  examples  shown  in  Figures  2  and  3*  Within  each  period 
parallelogram,  there  is  a  single  wave  peak.  It  occurs  in  a  region  of 
interaction  of  two  intersecting  "wave  crests".  The  crests  connect  adjacent 
peaks,  and  may  appear  as  ridges  of  constant  wave  amplitude,  as  in  Figure  2,  or 
merely  as  lines  of  steepest  descent  and  ascent  between  neighboring  peaks,  as 
in  Figure  3*  Ordinarily,  each  wave  crest  experiences  a  phase  shift  in 
crossing  a  region  of  interaction. 


39 


CONJECTURE  A:  Out  of  all  possible  KP  solutions,  the  only  ones  that  are 
(!)  real-valued,  (ii)  genuinely  two-dimensional,  and  (iii)  stationary  in  some 
uniformly  translating  coordinate  system  are  those  in  the  form  (1.9),  with 
genus  2. 

The  first  step  in  comparing  water  wave  data  to  KP  solutions  is  to 
transform  the  physical  data  into  KP  data  in  the  form  f(x,y,t).  This  is 
accomplished  via  (2.4),  with  F  -  0  in  (2.4b),  once  we  have  decided  on: 
(i)  the  magnitude  of  e;  and  (ii)  the  direction  of  the  x~axis  in  the 
horizontal  plane.  A  reasonable  definition  of  e  is  given  below  (2.3)*  A  more 
refined  definition  is  unnecessary,  because  adjustments  of  t  can  be  absorbed 
into  B  in  (2.8). 

The  direction  of  the  x~axis  should  be  approximately  the  principal 
direction  of  wave  propagation.  Small  readjustments  of  this  direction  can  be 
absorbed  into  a  in  (2.7),  but  one  can  expect  significant  errors  if  (2.5)  is 
used  to  model  water  waves  propagating  in  a  direction  much  different  from  the 
X-dlrection.  In  the  remainder  of  this  section,  we  assume  that  these  choices 
already  have  been  made,  and  that  the  water  wave  data  already  have  been 
transformed  via  (2.4)  into  KP  data. 

5.1  Wavelength  and  velocity  data 

As  discussed  in  §  3.  a  cnoidal  wave  is  determined  to  within  a  translation 
by  specifying  its  wavelength  and  speed  of  propagation.  The  KP  solutions  of 
genus  2  that  are  of  interest  here  are  two-dimensional  generalizations  of 
cnoidal  waves,  so  one  expects  that  a  similar  procedure  might  succeed  here  as 
well. 

Here  is  a  precise  statement  of  the  mathematical  problem.  Let  f(x,y,t) 
represent  a  particular  two-dimensional  KP  solution  of  genus  2.  Then  f  has  the 
form  (4.8),  and  is  specified  completely  by  the  eight  free  parameters  in  (4.1) 
and  (4.6).  Identify  eight  measurements  of  f(x,y,t)  which  are  sufficient  to 
determine  these  parameters,  and  give  an  explicit  algorithm  to  find  the 
parameters  from  the  measurements. 


The  practical  problem  differs  from  this  in  two  respects. 

(i)  Instead  of  f(x,y,t),  a  two-dimensional  wave  in  shallow  water  is 
observed.  The  wave  pattern  has  the  four  properties  listed  above  that  are 
common  to  all  two-dimensional  KP  solutions  of  genus  2,  and  we  use  Conjecture  A 
to  assert  that  the  water  wave  can  be  represented  by  a  two-dimensional  KP 
solution  of  genus  2. 

(ii)  Me  do  not  necessarily  require  that  the  eight  free  parameters  be 
determined  from  only  eight  measurements.  In  practice  it  is  often  convenient 

to  take  one  or  two  extra  measurements  and  to  minimize  the  effect  of 

measurement  errors. 

In  any  case,  how  does  one  find  experimentally  the  free  parameters  of  a 
two-dimensional  KP  solution  of  genus  2?  The  method  proposed  here  has  eight 
fundamental  steps. 

(1)  Find  a  period  parallelogram  at  a  fixed  time,  t. 

This  step  is  based  on  a  procedure  due  to  Arnol'd  (1978,  p.276).  Start 
with  any  point,  (x,y,t).  (From  here  to  step(5b  t  remains  fixed  at  t,  and  we 
shall  suppress  it.)  Because  f(x,y)  is  a  two-dimensional,  bi-periodic  function, 
there  is  an  infinite  array  of  points  in  the  plane  at  which  f  is  periodic  to 

f(x.y’);  i.e.,  at  which  f  and  all  of  its  derivatives  match  those  at  (x,y). 

This  collection  of  points  is  called  a  period  lattice  (relative  to  (x,y)). 
Pick  any  point  in  this  period  lattice  other  than  ("x, y") .  Denote  the  line 
through  these  two  points  by  ^  -  const.  Out  of  all  points  on  this  line  that 
are  in  the  period  lattice,  denote  one  closest  (but  not  equal)  to  (x,y)  by 

(*!  .y^. 

Out  of  all  points  off  of  this  line  that  are  in  the  period  lattice,  denote 
one  closest  to  the  line  by  (xg.yg)*  The  two  line  segments  connecting  (x,y)  to 
(x^y.,)  and  (x,y)  to  1x2^) i  are  two  sides  of  a  parallelogram  with  the 
following  properties. 

(i)  The  parallelogram  contains  exactly  four  points  of  the  period  lattice. 

These  four  points  lie  at  its  vertices.  (For  a  proof,  see  Arnol’d,  1978). 

(ii)  The  parallelogram  is  a  period  parallelogram  of  f;  i.e.,  f(x,y)  may  be 


defined  in  the  entire  plane  by  periodically  extending  the  function  defined  on 
the  parallelogram. 

(iii)  The  parallelogram  cannot  be  replaced  by  any  smaller  parallelogram 
without  losing  information  about  f(x,y). 

(iv)  The  enclosed  area  of  this  parallelogram  can  be  measured  directly,  and 
it  provides  information  about  the  parameters  in  question,  through  (4.16).  The 
method  described  here  does  not  use  this  area  explicitly,  but  it  can  be  used  as 
a  check  on  the  results  obtained. 

Define  ^(x.y)  as  in  (4.1),  and  choose  (p^,^)  so  that:  (1)  ^  -  0  along 
one  side  of  the  period  parallelogram;  (ii)  ^  -  2*  along  the  opposite  side. 
Define  in  the  same  way  using  the  other  two  sides  of  the  parallelogram. 

(2)  Find  the  wave  crests 

Within  each  period  parallelogram,  f(x,y)  has  one  complete  region  of 
interaction,  containing  the  maximum  of  f(x,y).  As  demonstrated  in  §4,  each 
region  of  interaction  is  connected  in  two  directions  (related  to  $3  and  4>i|  in 
(4.11))  by  saddles;  see  especially  Figures  3c,3d,3f»3g.  We  now  define  a 
"wave  crest"  to  be  a  straight  line  segment  along  a  saddle  connecting  two 
adjacent  regions  of  interaction.  These  are  wave  crests  in  the  sense  that 
along  any  other  line  that  intersects  one  of  these  line  segments  transversely, 
f(x,y)  attains  a  local  maximum  at  the  wave  crest.  In  this  paper,  "wave  crest" 
always  refers  to  a  saddle  between  two  peaks,  and  not  the  peak  itself. 

Emanating  from  every  region  of  interaction  are  two  pairs  of  parallel  wave 
crests.  We  will  eventually  associate  the  pair  possessing  the  larger  wave 
amplitudes  at  their  saddles  with  $3,  and  the  pair  possessing  the  smaller  wave 
amplitudes  with 

CONJECTURE  B:  The  wave  crests  in  each  period  parallelogram  are  defined  by 
$3  -  const,  and  -  const.,  where  $3  and  $4  are  defined  by  (4.11). 

(3)  Find  a  basic  period  parallelogram 
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The  period  parallelogram  found  in  step  (1)  is  not  unique.  Figure  5  shows 
two  dissimilar  period  parallelograms  on  the  same  period  lattice,  and  enclosing 
the  same  area.  Denote  the  sides  of  the  two  parallelograms  “  const, 

A  A 

and  by  ^ “  const,  respectively,  as  shown  in  Figure  5.  The  two  sets  of 
coordinates  in  Figure  5  are  related  by 

(l)  •  (: 

(5.1 ) 

as  the  reader  may  verify. 

More  generally,  let  (  ♦■) ,  $2 )  “  denote  the  sides  of  any  other 

period  parallelogram  on  this  same  period  lattice,  that  encloses  the  same  area 
as  that  of  ( ^ ^ )•  These  coordinates  are  related  by 

■  (:) 

(5.2) 

where  is  an  integer,  2x2,  unimodular  matrix  (i.e.,  integer  elements  and 
determinant  -  ±  1).  These  are  precisely  the  transformations  used  to  reduce  a 
real  Riemann  matrix  to  its  basic  form,  because  the  symplectic  transformations 
among  equivalent  Riemann  matrices  correspond  precisely  to  the  transformations 
among  period  parallelograms  with  the  same  enclosed  area  on  a  fixed  period 
lattice. 

The  two  kinds  of  transformations  relate  to  theta  functions  as  follows.  A 
real  2x2  Riemann  matrix,  plus  a  particular  period  parallelogram  on  a  lattice, 
determines  a  theta  function,  as  in  (H.3).  Suppose  the  Riemann  matrix  is 
transformed  to  an  equivalent  one  by  a  symplectic  transformation.  To  keep  the 
same  theta  function,  one  must  also  transform  the  period  parallelogram  by  the 
corresponding  integer  unimodular  matrix,  as  in  (5.2).  Thus,  for  a  given  theta 
function,  each  equivalent  Riemann  matrix  corresponds  to  its  own  period 
parallelogram  on  a  fixed  period  lattice.  (Here  and  in  what  follows  "period 
parallelogram"  refers  to  both  the  single  four-sided  figure  and  to  the  tiling 
of  the  entire  plane  by  this  figure.) 


Recall  that  a  real,  indecomposable,  2x2  Riemann  matrix  is  equivalent  to  a 
basic  Riemann  matrix,  for  which 

0<?n2ij  .  (5.3) 

Corresponding  to  the  basic  Riemann  matrix  is  a  particular  period 

parallelogram,  which  we  also  call  basic.  Here  is  how  to  identify  it.  Define 

♦•,(x,y)  and  $2(x,y)  as  in  (4.1),  so  that  ^f(x,y)/2x  takes  on  sequential 
integer  values  on  two  parallel  sides  of  the  parallelogram,  and  f2(x, y )/2lr 
takes  on  sequential  integer  values  on  the  other  two  sides.  The  basic  period 
parallelogram  is  identified  by  two  properties. 

(i)  It  is  related  to  the  parallelogram  in  step  (1)  by  an  Integer  unlmodular 
transformation.  In  particular,  it  encloses  the  same  area  as  the  parallelogram 
in  step  (1 ) . 

(ii)  The  wave  crests  are  related  to  the  lines  of  periodic!.:'  of  the  basic 
parallelogram  by 

♦3  (x.y)  -  <>,  (x,y)  -  I  ^  (x,y)  ,  (5.4) 

♦4  (x,y)  -  <t>2  (x,y)  -  A  +,  (x.y)  , 

where  (A,  A)  satisfy  (5.3)  and 

A  A  >  0  .  (5.5) 

He  show  in  Appendix  B  that  the  basic  period  parallelogram  is  uniquely 
defined  by  these  two  properties.  Here  are  some  other  properties,  which  follow 
from  these  and  from  the  results  in  14. 

(iii)  If  is  increased  by  2ir,  holding  <p2  fixed,  $3  changes  by  2 it. 
Conversely,  if  $2  is  increased  by  2it,  holding  $1  fixed,  $3  changes  by 
(-  2 u  A). 

(iv)  The  wave  crests  ($3)  associated  with  the  larger  amplitude  waves 


experience  a  smaller  phase  shift  (2«A),  while  the  smaller  amplitude  waves  feel 
a  larger  phase  shift  (2wA)  as  a  result  of  their  interaction. 

(v)  Given  a  period  parallelogram,  the  wave  crests  emanating  from  a  vertex 
are  constrained  by  (5.3)  to  lie  in  certain  two-dimensional  cones. 

(vi)  If  +2  lies  to  the  left  of  then  ^  lies  to  the  right  of  #2,  because 
of  (5.5).  The  interaction  does  hot  represent  a  common  rotation  of  both  wave 
crests. 

(4)  Measure  the  wave  numbers 

At  this  point,  we  have  identified  the  basic  period  parallelogram,  and 
have  located  the  two  sets  of  wave  crests  within  the  parallelogram.  Now  we  may 
measure  the  wave  numbers  (iij , w2»vi  »v2)  dir*ctly.  We  also  measure  A  directly. 

(i)  “  const  on  the  wave  crests  with  the  larger  amplitude  waves,  and  with 

the  smaller  deviation  from  the  nearest  line  of  periodicity. 

(11)  -  const  on  the  two  sides  of  the  parallelogram  most  nearly  aligned 

with  $3  -  const. 

(iii)  The  straight  lines  along  these  two  sides  are  defined  by 

UjX  ♦  v,y  ♦  Tiq  -  o,  p,x  ♦  v,y  ♦  J10  -  2ir  ,  (5.6) 

so  that  ( )  can  be  found  by  direct  measurement. 

(lv)  From  the  other  two  sides  of  the  parallelogram,  one  also  finds  (v2,v2) 
by  direct  measurement.  The  area  of  the  parallelogram,  of  course,  then 
satisfies  (4.16). 

(  v)  A  can  be  obtained  either  by  measuring  an  angle  and  using  (5.4),  or  by 
measuring  the  phase  shift  of  the  ^-crests  and  using  (4.25). 

(vi)  One  obtains  A  in  a  similar  way.  From  (4.11b), 

d/b  -  A/A  -  A2  •  (5.7) 

We  use  this  information  simply  as  a  check  on  (d/b),  but  it  is  not  necessary  to 


(5)  Measure  the  wave  speeds 


Because  the  wave  pattern  Is  a  two-dimensional  KP  solution  of  genus  2,  it 
is  stationary  in  a  uniformly  translating  coordinate  system.  Measure  this 
speed  of  translation.  Then  (4.17)  determines  (u^,u2). 

(6)  Deduce  the  Riemann  matrix 

Given  any  basic  Riemann  matrix  (b,d,A),  and  any  real  (114/^)  satisfying 
(4.13d),  then  (4.13a,b)  determine  ¥^,¥4.  Suppose  (¥j ,¥4)  were  generated  in 
this  way  from  some  (b,d,X,( V14/W1 ) ) . 

CONJECTURE  C:  Given  (1,  (ujj/pj),  ¥^,¥4),  (4.13a,b)  can  be  inverted  to  find 
(b,d)  uniquely. 

He  have  not  proven  this  assertion,  but  we  offer  the  following  evidence  in 
support  of  it. 

(1)  The  validity  of  the  corresponding  statement  for  genus  1  can  be  proven, 
(il)  He  have  tested  this  conjecture  numerically  by  first  generating  (¥^,¥4) 
from  a  particular  Riemann  matrix,  then  using  (4.13a,b)  in  a  simple  Newton-type 
root-finder  to  solve  for  (b,d),  given  (¥j  ,¥2* A, 114^1 ) .  Our  experience  has 
been  that  the  algorithm  always  converges  quickly  to  the  correct  values,  even 
if  our  starting  values  were  chosen  poorly. 

Returning  now  to  our  two-dimensional  KP  solution  of  genus  2,  in  effect  we 
have  already  measured  (Y^  ,¥4,  A,  (114/11^ )).  Then  it  follows  from  Conjecture  C 
that  (b,d)  are  determined  from  (4.13a,b).  Reasonably  accurate  starting  values 
can  be  obtained  from  (4.31),  (4.34)  and  (4.35). 

Conjecture  C  is  closely  related  to  another  conjecture  that  we  never  use 
explicitly,  but  which  is  implicit  in  this  entire  algorithm. 

CONJECTURE  D;  Every  bounded  real-valued  KP  solution  of  the  form  (1.9)  is 
generated  by  a  real-valued  Riemann  matrix. 


(7)  Check  Consistency 

(i)  The  ratio  (b/d)  has  been  determined  independently  in  steps  (4)  and  (6). 
These  must  agree. 

(ii)  Both  sides  of  (4.13c)  have  been  determined  independently.  These  must 
also  agree. 

(8)  Measure  the  phase  constants 

If  desired,  ^  r ,  $20  ln  can  be  measured  for  a  particular  coordinate 

system.  These  have  no  dynamic  significance.  In  any  case,  this  completes  the 
algorithm. 


5.2  Initial  Data 

A  second  way  to  specify  these  waves  is  with  initial  data:  to  specify  the 
elevation  of  the  water  surface  everywhere  at  some  fixed  time  (which  we  call 
t  -  0).  Suppose  f(x,y,0)  is  genuinely  two-dimensional,  is  periodic  in  two 
Independent  horizontal  directions,  and  is  given  completely  in  one  period 
parallelogram.  We  must  answer  two  questions:  (i)  Is  this  the  initial  data 
of  a  KP  solution  of  genus  2?  (li)  If  so,  how  does  the  wave  evolve  in  time? 
This  method  of  solution  was  first  given  by  Segur,  Finkel  &  Philander  (1983). 

Is  this  the  initial  data  of  a  two-dimensional  KP  solution  of  genus  2? 
Conjecture  A  is  not  effective  for  initial  data,  because  we  have  no  information 
about  the  time  evolution  of  the  solution.  Instead,  we  assume  that  f(x,y,0) 
has  the  form  (4.8)  and  derive  necessary  conditions  that  f  must  satisfy  in 
order  to  be  initial  data  of  a  two-dimensional  KP  solution  of  genus  2. 

For  any  function  satisfying  (4.8)  and  (4.1), 

■  U}  3^  +  U2&2’  «*y  "  ^i  ®i  +  v232,  (5.8) 

3^  "  3^  *  ui232, 


where 


47 


Now  integrate  f(x,y,0/  over  a  period  parallelogram  (with  area  A),  using 
(4.16): 

7  ( **1  v2  "  w2  V1 )  jj  dA  *  “  (5.9) 

ft  <*1  ft  d*2  (h31  +  “232  )  (3xtn0)  -° 

This  requirement  is  closely  related  to  the  normalization  in  (2.3). 

A  second  constraint  on  the  initial  data  is  obtained  by  multiplying  (1.8) 
by  (3xf),  and  integrating  over  a  period  parallelogram.  Most  of  the  integrated 


terms  vanish  identically. 

and  we  are 

left  with 

J  J dA*  | 

(3xf{x,y,0) 

)3  .o. 

(5.10) 

Similarly,  multiplying 

parallelogram  yields 

(1.8)  by 

3yf  and  integrating  over 

a  period 

//“*•  ( 

3xf(x,y,0) 

)  2  •  ^3yf (x,y,0)j  -  0. 

(5.11) 

It  is  clear  that  an  infinite  sequence  of  necessary  constraints  on  the  initial 
data  can  be  obtained  in  this  way.  It  is  also  clear  that  no  finite  number  of 
them  will  be  sufficient  to  guarantee  that  f(x,y,0)  is  the  initial  data  of  a  KP 
solution  of  genus  2.  Conditions  which  are  both  necessary  and  sufficient  are 
not  known  at  this  time. 

Consider  now  the  second  question.  Suppose  f(x,y,0)  is  two-dimensional,  is 
periodic  in  two  independent  directions,  and  is  given  on  a  period 
parallelogram.  Further,  suppose  that  f(x,y,0)  is  known  to  be  the  initial  data 
of  a  KP  solution  of  genus  2.  Find  f(x,y,t). 

Because  f(x,y,0)  defines  initial  data  for  a  two-dimensional  KP  solution 
of  genus  2,  the  initial  wave  pattern  simply  translates  with  some  uniform 
velocity  for  all  time.  Once  this  velocity  is  known,  the  initial  data  in  a 
period  parallelogram  plus  the  velocity  of  translation  determine  the  solution 
for  all  time.  There  is  no  need  to  find  the  underlying  Riemann  matrix  in  this 


case. 


The  velocity  of  translation  can  be  found  from  (4.17)  if  (  ,u>2  )  in  (4.1) 

are  known.  One  equation  for  (  u1  ,u)2  )  can  be  found  by  multiplying  (1.8)  by  f, 
and  multiplying  over  a  period  parallelogram.  The  result  is 

//«•  (Vt  *  6ffx  -  fx*  *  )  ■  «•  (5.12) 

It  follows  from  (5.1)  that 

( W1  v2  "  **2vl)  ft  *  (  “lv2  "  U2V1 )  fx  4  (**l“2  "  »*2«i)  fy  (5.13) 

Moreover,  ^  **i  ’  v2’ V1  ,v2  )  a11  can  be  measured  directly  from  the  given  initial 
data,  and  ii.^  *  u2vi  because  the  initial  data  is  two-dimensional. 
Substituting  (5.13)  into  (5.12)  yields 

//dA  (  Vx  *  »2fxry)  •  “1  *//«  (  »lfxfy  '  Vx  )  •  -j  (5.1«) 


♦ 


( 


W1V2  "  W2V1 


)  //d#  ( 


0. 


The  integrals  in  (5.14)  all  can  be  evaluated  at  t  -  0,  so  (5.14)  is  a/i 
algebraic  equation  for  ^  ,  u>2  )  • 

A  second  equation  for  (  , ui2  j  can  be  found  easily  if  (  j^^i  )  is 
rational,  which  we  now  assume.  Because  (  )  can  be  measured  from  the 

initial  data  only  to  a  finite  accuracy,  this  assumption  makes  no  practical 
limitations.  Its  consequence  is  that  f(x,y,0)  is  necessarily  strictly  periodic 
in  x  (holding  y  fixed).  Denote  this  x-period  by  L.  Because  f(x,y,0)  has  the 
form  (4.8)  by  assumption,  then  for  all  (  x0,y  ) , 


dx  f(x  +  xQ,  y,  0)  -  0. 


(5.15) 


Define 


i  (x,y;xQ)  -  if1  Jq  ds  •  sf  (x  ♦  xQ  ♦  a,y,0 )  (5.16) 

i  is  the  unique  anti-derivative  of  f  with  the  same  periodicity  as  f,  and  with 

zero  mean  in  xj  i.e.,  also  satisfies  (5.15)  and  f  -  3xtf.  The  corresponding 
anti-derivative  of  £  may  be  defined  in  a  similar  way. 

Now  multiply  (1.6)  by  that  anti-derivative  of  £  (i.e.,  the  second 

integral  of  f)  with  the  same  periodicity  as  £  and  f  and  integrate  over  a 
period  parallelogram.  The  final  result  is 

//»*  (  v2r2  -  v2f<y  )  •  ♦//«  (u,f<y  -  v,fS  )  •  »2  (5.17) 

♦  (<h'*  -  "a”! )  //“*  (3f3  -  'I  ♦  i  \  )  -  o. 

The  integrals  may  be  evaluated  at  t  •  0,  so  (5. 1*0  and  (5.17)  are  two 

linear  algebraic  equations  for  (  u1 ,u2  )  .  If  they  are  linearly  independent, 
their  common  solution  defines  (  ,u2  )  •  and  it  completes  the  specification  of 
the  given  KP  solution  of  genus  2.' 

CONJECTURE  E:  (5.1M)  and  (5.17)  are  linearly  independent,  and  they  define 
(<■>),  u2)  uniquely. 

It  is  easy  to  show  that  (5. 14)  and  (5.17)  are  linearly  independent  in  the 
small  amplitude  limit  (d  ■*  -  •,  b  ♦  -  •),  but  we  have  not  established  this 
property  in  general.  In  the  cases  we  have  tested  numerically,  (5.1 4)  and 
(5.17)  are  independent,  and  their  common  solution  agrees  with  that  in  (4.13). 
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Appendix  A 


Theorem ;  Let  ^  and  $2  be  independent  real  variables.  For  fixed  A  and  A 
(0  S  A2  2  A2  2  1/4),  define  +3  and  4>4  by  (4.11)  on  0  S  ^  <  2ir,  0  2  4>2  <  2*, 
and  by  periodic  extension  outside  this  rectangle.  Then  once  each  period 
63  -  const  and  -  const  experience  discontinuous  shifts,  whose  size  and 
direction  are  given  by  (4.25). 

Proof :  Figure  4  shows  a  period  lattice  in  ^  and  $2,  with  two  line  segments 

of  ^  -  const,  shifted  by  A^.  We  must  show  that  A$^  is  given  by  (4.25).  Let 
(x,y)  denote  Cartesian  variables,  and  let  ^  and  $2  be  related  to  (x,y)  by 
(4.1)  at  t  -  0.  Let  L  denote  the  normal  distance  between  $2  and  ($2  ♦  2ir). 
Then 


L  -  2 u22  ♦  v22)  “1/2 


(A.  1 ) 


Let  £1  denote  the  distance  along  ($2  -  const),  between  ^  and  (^  ♦  2ir),  as 
shown  in  Figure  4.  Their  product  gives  the  area  of  the  period  parallelogram: 


£,  L  ■  A 


, 

I  “1  v2  ‘  W2V1  I 


(A. 2) 


where  we  have  used  (4.16).  This  determines  £^ . 

The  unit  normal  vector  for  any  of  the  phase  lines 
♦j  -  const,  J  -  1,2, 3, 4  has  components  given  by 


ft*  - 


(A.3) 


'2 


The  angle  (Y)  between  the  two  lines,  <p2  m  const  and  ^  -  const,  is  given  by 


cos  Y  -  fU 


^4  + 


(y|  ♦  4 


y 


/2 


(4* 


v?)’ 


/2 


so 


sin  2  Y  -  1-  cos  2  Y 


(hv2  ~  ^2vl)2 
v22)  (u„2  ♦  v42) 


(A.4) 


Now  the  normal  distance  (ALjj)  is  one  side  of  a  right  triangle  whose  hypotenuse 
is  : 


ALjj  -  i1  sin  Y  - 

But  the  magnitude  of  the  wave-vector, 

ALjj  corresponds  to  a  phase  shift  of 

A$ty  ■  2irl  .  (A. 5) 

This  is  the  desired  result.  The  other  half  of  (4.25)  is  proved  by  similar 
means. 


_ 2itX _ 

2  2  \  1/2  * 

^4  +  v4  ) 

+  /  2  2  \  1  /2 

kjj  is  f  Mu  +  v]|j  1  ,  so  the  distance 
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Appendix  B 


Suppose  the  basic  period  parallelogram  of  a  fixed  wave  system  is 
identified  by  ^  ■  const,  and  $2  “  const.  The  wave  crests  are  identified  by 
♦3  ■  const,  and '  ^  -  const.,  where  ♦j  and  tj,  are  defined  by  (5. 4),  with  X  and 
X  satisfying  (5.3)  and  (5.5).  We  want  to  show  the  basic  period  parallelogram 
is  the  only  period  parallelogram  in  which  the  (fixed)  wave  crests  are  related 
to  the  lines  of  periodicity  through  relations  of  the  form  (5.3)-(5.5). 

Let  $a  ■  const,  and  $b  *  const,  along  the  sides  of  another  period 
parallelogram.  If  the  two  parallelograms  enclose  the  same  area,  then 
♦a,  <frb  are  related  to  4., ,  $2  by  an  integer,  unimodular  transformation  of  the 
form 


(B.1 ) 


In  this  Appendix,  let  on  ±  1 ,  with  om  and  on  independent.  Let  (a,b,c,d)  be 
integers  with 


Substituting  (B.1)  into  (5.4)  yields 
♦3 
*4 


a- Ac 

b-Ad 

c-Aa 

d-Ab 

(B.3) 


There  are  two  choices: 

Ci)  a  -  Ac  -  o2,  d  -  Ab  -  03,  b  -  Ad  -  A,  c  -  Aa  -  A;  or 


(ii)  c  -  Aa  -  o2,  b  -  Ad  -  03,  d  -  Ab  -  A,  a  -  Ac  -  A. 

For  either  choice. 


AA  >  0,  0<?iA2^, 


(B.4) 


We  will  show  that  in  either  case  one  is  forced  back  to  the  basic  period 
parallelogram. 

Assume  (i);  the  analysis  for  (ii)  is  similar,  and  we  will  not  repeat  it. 
Because  A A  >  0, 

(B.5) 


(B.6) 


a- Oj 
A 


d  -  0- 


b  ■ 


Therefore 


so  that 


A  -  c 


-  Aa  -  a|(  1-AA  j  -  o2J^  A  . 


o2  ♦  A  A 
1  -  AA 
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-VJ.  - 


.  A  ♦* 


It  follows  from  these  relations  and  (B.6)  that 

1  i  o2  a  S  |  .  (B.7) 

But  a  is  an  integer,  so  from  (B.7)  and  (B.5), 

a  -  02.  c  -  0.  (B.8) 

Similarly,  d  -  og,  b  »  0.  Therefore,  the  transformation  in  (B.1)  becomes 


(B.9) 

Any  choice  of  ( <>2*  °3 )  leaves  the  basic  period  parallelogram  intact.  This 
completes  the  proof  of  uniqueness  of  the  basic  period  parallelogram. 
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FIGURE  CAPTIONS 


Figure  1  .  Cnoldal  wave  solution  of  KP.  In  (3.2),  b  -  -  3.}  in  (3.1), 
P  -  0.5,  v  -  -  0.215,  »  -  -  0.325. 

Figure  2a.  Two  soliton  solution  of  KP.  In  (4.23),  hj  ■  x  ♦  My  -  M9t, 
ity  *  x  -  My  -  M9t,  exp  (A)  »  16/15. 

Figure  2b.  Genus  2  solution  of  KP  in  soliton  regime.  In  (M.6),  b  -  -  1, 
bX2  ♦  d  -  -  1,  X  -  0.15;  in  (M.1)  p,  -  p2  -  0.25,  v,  -  -  0.1263,  v2  -  0, 
u>i  •  -  0.9681  ,  «i>2  “  *“  0.2695. 

Figure  2c.  Contour  plot  of  wave  shown  in  Figure  2b.  This  plot  shows  clearly 
how  each  wave  crest  experiences  a  phase  shift  from  every  interaction  with 
another  wave. 

Figure  3  •  Genus  2  solutions  of  KP,  showing  some  of  the  variety  of  wave  forms 
available. 

Figure  3a.  One  wave  is  dominant,  b  -  -  2,  d  -  -  3,  X  -  0.3,  Hi  -  p2  ■  0.5, 
Vj  ■  -  0.316M,  v2  ■  0,  <■>]  -  -  1.227,  w 2  m  0.063. 

Figure  3b.  Both  wave  crests  are  evident,  b  •  -2.5,  d  ■  -  3«,  X  -  0.3, 
Pj  ■  P2  *  0.5,  v1  -  -  0.2579,  v2  ■  0,  -  -  0.7360,  ii>2  *  0.1313. 

Figure  3c.  A  comparison  of  the  waves  here  and  in  Figure  3b  shows  the  effect 
of  the  scaling  symmetry  in  (2.8).  Here  b  -  -  2.5,  d  »  -  3»»  X  -  0.3, 
W1  •  p2  ■  0.8,  V|  ■  -  1.056,  v2  -  0,  -  5.56M,  «£  *  **  0.55M9. 

Figure  3d.  Contour  plot  of  the  wave  shown  in  Figure  3c.  Here  a  basic  period 
parallelogram  is  shown,  along  with  the  wave  crests  corresponding  to  $j( —  ) 
and  ^(—  •— ).  Note  that  the  wave  crests  do  not  lie  along  the  directions  of 
periodicity. 


Figure  3e.  Contour  plot  of  wave  pattern  in  the  small  amplitude  limit, 
b  ■  ~  8,  d  •  -  8,  A  ■  0.1,  U]  m  V2  m  ^  •  vl“”  1*808,  V2  ■  0,  <*)■)■“  8.820, 
«2  •  0.9837.  In  this  limit  (and  only  here)  the  wave  crests  align  with  the 
directions  of  periodicity. 

Figure  3f.  "Typical"  genus  2  solution  of  KP,  away  from  any  limiting  case, 
b  -  -  d  -  -  3.8,  A  ■  0.5,  Ui  ■  1*2  “  ^ '  V1  “  “  0.7147,  V2  ■  0,  uj  ■  -  1.300, 
U2  **  0.0885. 

Figure  3g>  Contour  plot  of  wave  pattern  in  Figure  3f.  The  wave  crests  and 
phase  shifts  are  marked. 

Figure  4.  Period  diagram  of  hypothetical  wave  pattern,  showing  lines  of 
periodicity  (tj  and  *  wave  crests  and  phase  shift.  The  same  effects  in  an 
actual  wave  pattern  are  shown  in  Figure  3g. 

Figure  5  .  A  fixed  period  lattice  supports  many  period  parallelograms  that 
enclose  the  same  area.  Two  are  shown  here. 


